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to  setting  up  a  systematic  approach  for  deriving  state  and  measurement  models.  Properties  of 
quadratic  and  bilinear  vectors  are  investigated  and  relationships  between  these  and  linear 
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Executive  Summary 

A  number  of  important  control  and  estimation  problems  in  the  field  of  aerospace  and 
avionics  involve  dynamical  models  that  are  quadratic  and  bilinear  functions  of  system 
states  and  inputs.  Examples  of  some  of  these  problems  are;  target  tracking,  vehicle 
navigation,  guidance  and  autopilot  design.  While  there  are  a  number  of  different  ways 
of  posing  (modeling)  the  above  problems,  it  appears  that  by  appropriate  choice  of  state 
and  input  variables,  systems  involving  translational  and  rotational  kinematics  can  be 
modeled  as  a  set  of  differential  equations  that  contain  linear  plus  quadratic  and/or 
bilinear  expressions.  The  well-known  Riccati  differential  equation  that  occurs  in  the 
synthesis  of  optimal  control  and  estimation  problems  is  also  of  this  class  of  differential 
equations. 

In  this  report,  formal  definitions  of  free  and  forced  quadratic/ bilinear  dynamical 
systems  are  given  and  the  algebraic  structure  of  this  class  of  systems  is  explored  with  a 
view  to  setting  up  a  systematic  approach  for  deriving  state  and  measurement  models. 
This  class  of  dynamical  systems  is  characterized  by  a  set  of  first  order  differential 
equations  with  the  RHS  that  contains  linear  and  quadratic  terms  in  system  states  as 
well  as  bilinear  terms  involving  state  and  input  (or  control)  variables.  Properties  of 
quadratic  and  bilinear  vectors  are  investigated  and  relationships  between  these  and 
linear  vectors  are  established.  Systematic  procedure  for  constructing  quadratic  state 
and  biUnear  state-input  vectors  is  derived.  The  quadratic  algebraic  structure  is  applied 
in  the  formulation  of  a  state-space  model  for  a  second  order  approximation  of  a  general 
(analytic)  non-linear  system. 

The  concept  of  quadratic  and  bilinear  generator  matrices  is  introduced  that  allows 
linear  (state  and  state-input)  vectors  to  be  mapped  onto  quadratic  and  bilinear  vectors. 
Properties  of  these  generator  matrices  are  explored  and  it  is  shown  that  corresponding 
inverse  generator  matrices  may  be  defined  that  allow  quadratic  and  bilinear  vectors  to 
be  mapped  onto  linear  state  and  linear  input  vectors.  State-space  representation  of  the 
system  output  or  measurement  model  is  also  derived  for  the  case  where  this  contains 
linear  and  quadratic  terms  in  states  and  bilinear  terms  in  state-disturbance  (noise) 
variables. 

The  state-space  structure  of  quadratic/biLinear  dynamical  systems  considered  in  this 
report  should  facilitate  analysis  of  this  class  of  non-hnear  systems  and  possibly  lead  to 
general  synthesis  techniques  or  extension  of  linearised  techniques  to  these  problems. 
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1.  Introduction 

A  number  of  important  control  and  estimation  problems  in  the  field  of  aerospace  and 
avionics  involve  dynamical  models  that  are  quadratic  and  bilinear  functions  of  system 
states  and  inputs.  Examples  of  some  of  these  problems  are:  target  tracking,  vehicle 
navigation  and  guidance  and  autopilot  design  [1-8].  While  there  are  a  number  of 
different  ways  of  posing  (modeling)  the  above  problems,  it  appears  that  by  appropriate 
choice  of  state  and  input  variables,  systems  involving  tanslational  and  rotational 
kinematics  can  be  modeled  as  a  set  of  differential  equations  that  contain  linear  plus 
quadratic  and/or  bilinear  expressions.  The  well-known  Riccati  differential  equation 
[9,10]  that  occurs  in  the  synthesis  of  optimal  control  and  estimation  problems  is  also  of 
this  class  of  differential  equations. 

Other  applications  of  bilinear  systems  are  in  die  field  of  industrial  process  control, 
economics,  biological  systems,  and  chemical  process  modelling  and  control  [11-13].  In 
addition  to  the  applications,  various  researchers  have  also  reported  on  important 
developments  in  bilinear  systems  theory  and  related  stabilization,  observability  and 
controllability  issues,  including  the  Lie-algebra  approach  [14-16].  The  control  synthesis 
problem,  involving  bilinear  systems,  has  been  addresses  by  a  number  of  authors 
utilising  the  optimum  control  theory  [17-20].  While  the  above  research  lays  down  a 
fundamental  and,  in  many  respects,  rigorous  basis,  for  system  analysis  and  synthesis 
for  this  class  of  non-linear  systems,  there  still  remains  a  need  to  apply  the  results  of 
these  investigations  to  high  dimension  systems  (involving  large  number  of  state  and 
control  variables)  that  arise  in  many  practical  engineering  problems.  Motivated  by  this 
need,  the  current  report  addresses  issues  relating  to  the  algebraic  structure  of  quadratic 
and  bilinear  dynamical  systems  and  proposes  an  approach  for  state-space  formulation 
of  the  dynamical  models  that  should  allow  for  systematic  and  generalized  procedures 
in  dealing  with  practical  control  analysis  and  synthesis  problems. 

In  this  report,  formal  definitions  of  free  and  forced  quadratic/ bilinear  dynamical 
systems  are  given  and  the  algebraic  structure  of  this  class  of  systems  is  explored  with  a 
view  to  setting  up  a  systematic  approach  for  deriving  state  dynamics  and  measurement 
models.  This  class  of  dynamical  systems  is  characterized  by  a  set  of  first  order 
differential  equations  wiih  the  RHS  that  contains  linear  and  quadratic  terms  in  system 
states  as  well  as  bilinear  terms  involving  state  and  input  (or  control)  variables. 
Properties  of  quadratic  and  bilinear  vectors  are  investigated  and  relationships  between 
these  and  linear  vectors  are  established.  Systematic  procedure  for  constructing 
quadratic  state  and  bilinear  state-input  vectors  is  derived.  The  quadratic  algebraic 
structure,  developed  in  this  report,  is  applied  in  the  formulation  of  a  state-space  model 
for  a  second  order  approximation  to  a  general  (analytic)  non-linear  system. 
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The  concept  of  quadratic  and  bilinear  generator  mati  ices  is  introduced  that  allows 
linear  (state  and  state-input)  vectors  to  be  mapped  onto  quadratic  and  bilinear  vectors. 
Properties  of  these  generator  matrices  are  explored  and  it  is  shown  that  corresponding 
inverse  generator  matrices  may  be  defined  that  allow  quadratic  and  bilinear  vectors  to 
be  mapped  onto  linear  state  and  linear  input  vectors.  State-space  representation  of  the 
system  output  or  measurement  model  is  also  derived  for  the  case  where  this  contains 
linear  and  quadratic  terms  in  states  and  bilinear  terms  in  state-disturbance  (noise) 
variables. 

The  state-space  structure  of  the  quadratic/bilinear  d)mamical  systems  considered  in 
this  report  should  facilitate  analysis  of  this  class  of  non-linear  systems  and  possibly 
lead  to  general  synthesis  techniques  or  extension  of  linearised  techniques  for 
applications  to  these  problems. 

The  algebraic  structure  of  free  quadratic  d5mamical  systems  is  considered  in  section  2 
of  this  report  and  includes  consideration  of  various  properties  of  quadratic  state 
vectors  and  their  relationships  with  linear  state  vectors.  Bilinear  vectors  are  considered 
in  section  3  along  with  the  properties  of  these  vectors  and  their  relationships  with 
linear  vectors.  Sections  4  and  5  consider  forced  dynamical  systems  and  measurement 
models  where  the  latter  includes  linear  as  well  as  quadratic  terms  in  state  variables  and 
bilinear  terms  involving  state  and  disturbance  variables.  In  section  6  general  non¬ 
linear  dynamical  systems  with  analytic  non-linearity  is  considered  and  a  second  order 
perturbation  state-space  model  of  the  system  is  derived  in  a  quadratic/bdinear  vector 
form. 


2.  Free  Dynamical  System 

In  this  section  we  formalise  the  definition  of  a  class  of  non-linear  dynamical  systems 
whose  time  evolution  of  states  may  be  expressed  via  a  set  of  first  order  differential 
equations,  the  right  hand  side  (RHS)  of  which  contain  botii  linear  and  quadratic  state 
terms.  A  formal  definition  will  lead  to  a  state  vector  characterisation  (or  a  state-space 
representation)  of  these  systems.  The  class  of  dynamical  systems  being  considered  in 
this  section  will  be  referred  to  as  free  (or  uncontrolled)  quadratic  dynamical  systems. 

2.1.  Definition 

An  n-state  non-linear  dynamical  system  will  be  defined  as  a  “free  quadratic  dynamical 
system"  (free-QDS)  if  its  dynamics  can  be  expressed  by  the  following  set  of  differential 
equations: 


dt 


Xiit)  =  X  'Zdj  i  =  . n 


(2-1) 


7=7  k=Ij=k 
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Where: 

— (•) :  is  the  time-derivative 
dt 

{x  j {t);  j  =  are  linear  system  states 

{xjt  {t)xj  (i);  k  =  1,2,3,...,  n;  j  =  k  +  l,k  +  2,...,n\:  are  quadratic  system  states 

{fl,-  f  (i);  i  =  l,2,3,...,n;  j  =  l,2,3,...,n\:  are  linear  system  state  coefficients 

^k'j  "  l,2,3,...,n;  k  =  l,2,3,...,n;  j  =  k,k  +  w}:  are  quadratic  system  state 

coefficients 

Alternatively,  the  free-QDS  of  equation  (2-1)  may  also  be  written  as: 

W  +  (0;  j  =  I2,...,n  (2-2) 

j-j 

Where: 

l  =  +  ^  =  {k- ^ j  +  7 

i  if  k  =  j 

^k.j  = 

yi  if  ki^  j 


q['j{t)=akjbij{t) 


Remarks: 


1.  The  RHS  of  equations  (2-1)  and  (2-2)  contain  linear  terms  {x,  (f)}  as  well  as  quadratic 
terms  jx^^  (i)x^  (i)}.  The  coefficients  of  the  quadratic  terms  }  are  related  to  the 

elements  of  a  symmetric  matrix  that  may  be  used  to  generate  the  scalar 


quadratic  function  |x^  x 

x  =  [x;  X2  x„]^[l,2].  Here: 


from  a  "linear  state  vector"  (Isv): 


if  k  =  j 
if  k^j 
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2.  The  RHS  of  equation  (2-2)  is  identical  to  that  of  equation  (2-1)  except  for  the 

coefficients  of  the  quadratic  terms.  The  terms  have  been  replaced  by  the 

terms  fbj,  (r)}  in  equation  (2-2).  This  substitution  has  been  made  to  allow  for 
state  differential  equations  to  be  expressed  in  matrix  notation,  as  will  be  seen  later. 
The  reason  for  inti  oducing  the  normalisation  factor  j  =  V2  or  l\  is  explained 
later  in  this  section. 

3.  In  expanded  form  the  differential  equation  (2-2)  may  be  written  as: 


4- X,. (/)  =  o,- i{t)xi [t) + fl,. , {t)x2 (0 + - ■  {t)x„ (/)  +  bjj {t)xj- (t) + ^l2bfj {t)xj {t)x2 {t) 

at 

+V2Z>,.  3  (i)  +. . . .  +  VJZ),-  „  {t)x ,  {t)x„  (/) + 6,  (i)A- (?) + ^^2bj„+2  (0^2  (0^5  (0 

W2Z),.  {t)x2  {t)x^  (/)+....  +  {t)x2  {t)x„  (i)  +  ....  +  6  r  j  W 

’  2 

for  {i  =  J,2,....,n}  (2-3) 


+V2/>  r  «(„+]) 


i{t)x„_jit)x„{t)+b 

■«(«+?)■ 

1, 

2 

4.  Note  that,  in  the  quadratic  part  of  the  above  equations,  there  are 
coefficients  {&, /(?)},  for  each  i,  associated  with  the  quadratic 
terms  {a^.  {t)xj  (i)}. 


5.  The  introduction  of  the  factor  Vi  in  the  coefficients  was  suggested  in  [3],  and  as 
will  be  shown  later,  it  lends  a  useful  property  to  the  "metric"  for  a  "quadratic  state 
vector"  (qsv).  However,  before  we  discuss  this,  it  will  be  instructive  to  consider  an 
example  of  a  second  order  free-QDS.  The  argument  (?)  is  suppressed  in  the  sequel, 
it  being  understood  that  variables  and  coefficients  are  functions  of  time  unless 
otherwise  stated. 


2.2.  Example 

This  example  has  been  selected  in  order  explore  some  useful  properties  of  a  qsv  and 
relate  these  to  the  more  familiar  properties  of  an  Isv.  We  consider  a  free-QDS  for  n=2; 
in  this  case  the  system  differential  equations  may  be  written  as: 
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— X]  =  Qj  ]Xj  +  Oj  iXt  +  bj  jxf  +  42b,  2^1X2  +  b, 
dt  ' 

—  X2  =  OtjX]  +  02  2^2  +  b2jX,~  +  42b2  2X1X2  +  b2jX2~ 


Or,  in  matrix  notation,  differential  equation  (2-4)  may  written  as: 


^  ^1,2  2Cj  b,,  b]2  b,2  X] 


dt  X9J  \j22J  (2 2, 2  X]  b2i  b2  2  ^2,5 J  4~2X]X 


Equation  (2-5)  is  of  the  form: 

dt 

Where: 

x^^^  =  [xy  :  is  a  [2x7]  linear  state  vector  (Isv) 

x^~^  =  [xy"  42xjX2  X2~J  :  is  a  [2x7]  quadratic  state  vector  (qsv) 

[,4^^^]=  :  is  a  [2x2]  linear  state  coefficient  matrix 

1^2,1  ^2,2} 

r  1  ~  b  b  b  ~ 

[5/^^]=  bt  1.2  1,3  .  is  a  [2x2]  quadratic  state  coefficient  matrix 

^27  ^2,2  ^2,7 


Remarks: 

1.  For  this  example,  the  quadratic  state  vector  (qsv)  may  be  regarded  as  a 

vector  spanning  the  space  whose  elements  are  the  terms  of  a  homogeneous 
quadratic  polynomial  in  (xy,  x^}  arranged  in  lexicographic  order. 

2.  The  inclusion  of  the  factor -v/2  imparts  a  desirable  property  to  the  metric  of  the  qsv. 
Thus,  analogous  to  tire  definition  of  a  metric  for  an  Isv  let  us  define  a  metric  for  a 
qsv  as  the  square  root  of  fire  sum  of  squares  of  the  quadratic  elements.  That  is,  for 
n=2: 


II  =  Is  P I '  =  I  I”  +  [42x1X2  p  +  (x3^  )■  y  =  \c4  +  X2~  ]=  |U^^^  II  (2-7) 


3.  It  is  interesting  to  note  that  by  defining  the  qsv  as  we  have  done  above,  and  by 
adhering  to  the  familiar  definition  of  a  metric  in  a  linear  vector  space,  a  relationship 
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between  the  metrics  of  the  qsv  and  the  /si'  may  be  established.  That  is,  the  norm  of 
the  qsv  is  equal  to  the  square  of  the  norm  of  the  associated  Isv. 

4.  A  qsv  formed  by  including  the  factor  VJ  will  be  referred  to  as  a  'normalised'  qsv  or 
simply  as  a  qsv.  The  'non-normalised'  qsv,  that  is  =  [xy"  XyX?  J  can  also 
be  used  in  the  formulation  of  a  free-QDS.  However,  in  this  case: 


I 


5.  We  observe  that  a  second  order  quadratic  polynomial  of  the  form: 

Xy"  +  VJx;X2+X9^ 


(2-8) 


(2-9) 


whose  terms  form  elements  of  the  qsv  x  ^  ,  may  be  expressed  as  a  combinations  of 
scalar  quadratic  functions  (forms)  using  the  Isv  x^^^  as  follows: 


..JO)-  yo  j  1  ix'o  f  )+ y-i  j  [e,,  i  (,» ) 

[k=lj=k  J 


Where: 


J  0 
0  0 


1 

V2 


0  1 
1  0 


0  0 
0  1 


(2-10) 


(2-11) 


6.  We  have  selected  these  [^^yj  matrices  to  be  symmetric.  These  matrices  serve  as  the 

basis  for  the  quadratic  pol3momial;  in  fact,  any  second  order  quadratic  polynomial 
may  be  written  (in  terms  of  a  linear  combination  of  matrices  such  as  \_E^j\)  as 

follows: 


>;=‘‘scalay  (2-12) 

{k=Ij=k  J 

7.  The  basis  matrices  \e^j\  may  also  be  used  to  construct  the  qsv  from  the  Isv,  which 

allows  us  to  establish  a  relationship  between  the  qsv  and  the  associated  Isv  as 
follows: 
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1 

_ J 

i[Eu]x 

I.J 

II 

42XjX2 

= 

^\Ei,2\^ 

1 

_ 1 

i\E2,2]x 

Xj 

X-, 

0 


0 

IL 

r2 

X2 


Xj 

X2 


(2-13) 


Where: 


0 


X2  Xj 

V2 

0  X2 

is  a  [3x2]  matrix 


/[Eu] 

/fe.j] 


\Ejj\_  [Ej_2\c  [£,,]xf 


(2-14) 


8.  The  matrix  \,  which  will  be  referred  to  as  the  ''quadratic  generator  matrix" 
(qgrn),  maps  the  Isv  x^^^  onto  the  associated  qsv  x^'^ .  It  will  be  shown  later  that 
provided  x^^^  ^0,  there  exists  the  "inverse  quadratic  generator  matrix"  (iqgm): 

that  maps  the  qsv  x^~^  onto  the  associated  Isv  x^^^ .  (Note  the 

change  between  the  sub  and  super-scripts  between  qgm  and  iqgm  symbols).  This 
inverse  relationship  follows  from  equation  (2-13): 

Now,  x'-^  =  [xj]j]x^’^ 

i-e :  x^‘^  =[xl2j\x^~^  (2-15) 


Where: 


:  wiU  be  referred  to  as  an  inverse  quadratic 


generator  matrix  iqgm. 


In  fact,  using  equation  (2-14)  and  after  some  algebraic  manipulation  we  get: 


(2-16) 
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and  the  [2x5]  iqgm  is  given  by: 


[2]  \-  ^ 


^  W  {  f  2  \2ir  2 

X,'  7  X->  Xj  1  Xj  X-,  XjX7 

b?  hivf 


’ll/'  2  \  2 

Xj'X-,  \  Xi  -)  X-,  XjX-> 


9  ^9 

X-,  x,'  +-^ 
-  ^  2 


(2-17) 


With:  /I  = 

2 

Note:  that  A  ^  0;  provided  x^O. 


9.  The  results  of  this  example  are  generalized  in  the  sequel;  however  before  we  leave 
this  example,  we  shall  use  equation  (2-13)  or  the  qgin,  to  give  an  alternative  state 
space  representation  of  the  free-QDS  analogous  of  equation  (2-6).  That  is: 


(2-18) 


2.3.  State-Space  Representation  of  a  General  n-State  Free-QDS: 

In  vector  notation  a  general  n-state  free-QDS  [see  equation  (2-6),  Example  2.2]  may  be 
written  as: 


(2-19) 


Where: 

=  [x^  X;  ....  x„ ]^ :  is  an  [w x /]  linear  state  vector  {Isv) 

yf2XjX2  42xjX2  ....  42xjX„  Xo^  42X2X2  42x2X^  ... 

X3  '42Xfj_2Xjj  ^n-l"  42x„_]Xu  Xfj~ 

:  is  an  quadratic  state  vector  (qsv) 

Note  that  the  l‘<‘  term  {xf  }  of  the  vector  x^^^  is  given  by: 


xj^^  =aj^jX^Xj;  k  =  j  =  k,k  +  l,k  +  2,....,n 


l  =  {k-]in—-\  +  j;  akj  = 


1  ifk  =  j 


V2  if  kit  j 
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^11  ^12  ^1,2  "*■  ^1 


O-ij  O-y^y  •••* 

k"]= .  :  is  an  [wXw]  linear  state  coefficient  matrix 


^n.l  ^n.2  ^i>,3  ■"■  ^»,h_ 


^LI  ^1,2  ^1,3 

^2,1  ^2,2  ^2,3 

[«'-’']=  . 

b„.i  K2  K3 


)(»-/)  j  »(»+/)  j 

'1  J  '1  J  .  r  w(«+/)]  , 

:isan  nx—^ — ^  quadratic-state 


coefficient  matrix 


Utilising  the  concept  of  the  qgm  introduced  in  Example  (2-2),  we  may  write  an 
alternative  form  of  a  general  n-state  free-QDS  as: 


(2-20) 


Where: 


[a" (;/  J :  is  the  ^  x  n  qgm  defined  in  equation  (Bl-1),  Appendix-B. 


Remarks: 


1.  We  make  the  observation  that  provided  x^O,  none  of  the  columns  of  a  qgm  ] 
equation  (B2-1)  may  be  expressed  as  linear  combinations  of  the  other  (n-1)  columns; 
that  is,  it  has  a  rank=n,  and  hence  the  matrix  is  of  full  rank  and 


therefore  invertible. 


2.  The  qgm  for  a  general  n-state  system  may  be  constructed  using  the  basis  matrices 
introduced  earlier  [see  equations  (2-10)  to  (2-14)].  It  follows  from  equation  (2-13) 
that 


kwl=fc,/k  [^7,2 k  Kik  fc,«k  fc,2k  Kvk  fe.^k  •••• 
fe„k  fc,5k  •••■  fc-.v,k  k,-7,»-ik  k,-7.«k  fc.„kr 


(2-21) 


9 


DSTO-TR-1497 


Where: 


'  I 

0 

0  ... 

o' 

0 

1 

0 

...  o' 

'  0 

0 

0 

]' 

0 

0 

0  ... 

0 

1 

0 

0 

...  0 

f 

0 

0 

0 

0 

k;]= 

0 

0 

0  ... 

0 

;k..r 

1 

0 

0 

0 

...  0 

k_H  ]— 

7 

7^ 

0 

0 

0 

0 

0 

0 

0  ... 

0 

0 

0 

0 

...  0 

1 

0 

0 

0 

'o 

0 

0  ... 

o' 

'0 

0 

0 

...  o' 

0 

0 

0 

o' 

0 

1 

0  ... 

0 

0 

0 

1 

...  0 

0 

0 

0 

0 

kv’]= 

0 

0 

0  ... 

0 

/■  \^2,3 

I 

0 

1 

0 

...  0 

0 

0 

0 

0 

0 

0 

0  ... 

0 

0 

0 

0 

...  0 

0 

0 

0 

1 

Note  that  tlie  basis  matrices  \  have  elements  given  by: 


for  j  =  k 
for  j  It  k 


1  for  r  =  s 
0  otherwise 


—p=  for  r  =  k,s  =  j  and  r  =  j,s  =  k 

'  42 

0  otherwise 


2.4.  Properties  of  the  General  n-State  Quadratic  Generator  Matrix 
(qgm): 

It  can  easily  be  verified  that  for  the  n-state  qsv  and  the  corresponding  qgm,  the 
following  properties  hold  (the  proof  can  be  arrived  at  in  a  manner  similar  to  that  of  the 
2-state  case,  and  is  left  to  the  reader): 


2.4.1.  Metric  Property: 


J2J 


/«(«+/) 


I 

1=1 


2.4.2.  Mapping  of  Isvxf'^  onto  qsv  through  qgm  \^pj 


[2] 


k=I 

J=k 


(2-22) 


(2-23) 
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2.4.3,  Inverse  mapping  of  qsv  ontolsv  [X|fjjJ; 

(2-24) 

Where:  k,'i;)={  [jr'//f  ^/;/]}  '\xfff 

2.4.4.  The  Bilinear  Property: 

Given  any  two  "linear  vectors"  e  R'' ,  the  quadratic  generator  matrix  (qgm) 

k/f;  ^  defined  as  in  equation  (B2-1)  will  be  called  the  "bilinear  generator 
matrix"  (bgm)  when  it  operates  on  e  R"  to  yield  the  "bilinear  vector" (bv) 

n(n+l)^ 

Similarly,  the  quadratic  generator  matrix  (qgni)  pjfj  Je  i?  -  defined 
as  in  equation  (B2-1),  (with  x,  replaced  by  uf  will  also  be  called  the  "bilinear  generator 
matrix"  (bgm)  when  it  operates  on  7?”  to  yield  the  bilinear  vector 

Where: 

XjUj  -^{xjU2+X2ttj)  -^{xjU^+X^Uj)  ....  -^{xjU„  +X„Uj)  X2U2  {x  2U  j  +  X  2) 

42  \2  V2  V2 


-^{X2U„+X„U2)  .. 


(2-25) 


The  vector  (.,.)  as  defined  above  in  (2-25),  includes  terms  in  sum  of  products  of  the 
elements  {x4,m^,^  =  7,2, and  7  =  A:, A: +  7, taken  in  lexicographic  order.  This 
vector  will  be  referred  to  as  the  bilinear  vector  (bv)  of  tihe  corresponding  linear  vectors 
Note:  that  the  order  of  the  arguments  in  {.,)  is  unimportant;  i.e. 
are  interchangeable.  Also,  since  both  are  n-order  vectors,  we 


shall  classify  this  bv  as  having  the  dimension 


«(n  +  7) 


x7L  Clearly: 


(2-26) 


Remarks: 
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1.  In  view  of  tlie  properties  (2.4.2)  and  (2.4.3)  above,  we  shall  formally  recognize  the 
bgms  [^ffj  I  Ppj  Je  ^  '  as  transformation  matrices  that  define  the  mappings: 

d.  It  can  easily  be  verified  that  the  norm  of  the  bilinear  vector  satisfies  the 
relationship: 


7') 


n] 


I  + 


Jt] 


.m 


1  f,]"^ 

—  x'  ^ 


Where:  y  is  the  angular  separation  between  the  vectors  and  . 


2.  It  will  be  instructive  at  this  stage  to  consider  a  bilinear  vector  where  the  dimensions 
of  the  two  linear  vectors  involved  are  different.  Accordingly,  given  any  two  linear 

vectors  e  R” e  R”' ;m  <n ,  the  corresponding  bilinear  generator  matrices 
(bgms)  are  given  by: 


—[2] 
^  [1] 


R 


lx  m 


Jlxm 


u 


[2] 

fU 


j^txn.j  m{2n-m  +  }) 


-ilxn  2 

These  matrices  are  defined  in  equations  (Bl-2)  and  (Bl-3),  Appendix  B. 


This  result  along  with  the  others  given  in  this  section  will  be  used  for  developing 
the  concept  of  a  "bilinear  state-input  vector"  (bsiv)  and  of  forced  (or  controUed)- 
quadratic/bHinear  dynamical  systems  (Forced-QBDS). 

2.4A.1.  Example 


For  n=3,  m=2,  the  bgms  are  given  by  the  following  relationships: 
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Xj  0 


^2  42 
M,.,  -Ti  ■ 


0  X-, 


Uj  0  0 

J!2_  JiL  n 

42  42 

0  0 

42 

0  u->  0 


0  0^ 


XjU]  -^{xoVj  +XJU2)  -^{x^U])  X2U2 

42  42  V-? 


Remarks 

1.  The  "bilinear  n-state,  m-input  vector"  (bsiv)  is  given  by: 


XjUi  -^{xjUt +X',Uj)  -^{xiUj+X^Uj)  ...  -^{x„Uj)  XyUj 

42  '  ~  42  42 


{x 2II 2  +  X 2^  ...  ^j—ix^U2^  ...  ip^n^ni-l^ 


(2-27) 


(2-28) 


2.  Finally,  we  end  this  section  by  giving  a  block  diagram  for  a  free-QDS  as  reflected 
by  the  state  vector  equation  (2-19).  This  is  given  in  Figures  2.1. 

3.  For  sake  of  convenience  the  sub-scripts  outside  the  square  brackets  of  the  hgms 
have  been  dropped,  that  is: 


Figure  2.1  Block  Diagram  of  the  free-QDS  ofecpiation  (2-19) 
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3.  Bilinear  State-Input  Vector 

Before  we  consider  forced-QBDS,  we  shall  formalize  the  sti  ucture  and  definition  of  a 
"bilinear  state-input  vector"  {bsiv)  generated  from  two  linear  vectors,  namely,  the  n- 
dimensional  state  vector  e  R"  and  the  m-dimensional  input  vector  e  R'" .  In 
order  to  develop  this  concept  we  shall  utilize  the  familiar  scalar  bilinear  forms  and  its 
various  properties.  Accordingly,  we  consider  the  following  scalar  bilinear  form: 

77,4w,J'1)= 

k=lj=l  ■  k=J_H  ' 


In  a  more  familiar  matrix  notation  equation  (3-1)  may  be  written  as: 


k=lj=l 


Where: 

[5'^'^]=  .•  k  =  -  l,2,...,n;  with  m  <  «|:  is  an  [r  x  /??]  matrix 

;  j  -  7, 2, is  an  [«  x7]  Isv 

.•  j  =  7,2,...,  w}:  is  an  [w;  x  7]  "linear  input  (or  control)  vector"  (liv) 


(3-2) 


By  expanding  the  double  summation  in  equation  (3-2)  and  rearranging  the  terms,  it  is 
easily  verified  that: 


k=lj=l  k=IJ=k 


(3-3) 


Where: 


^k,J  ] 


for  j  =  k 


1  otherwise 


Since  Uj  =  0  for  j>m,we  may  w.Lo.g  assume  =  0  for  j>m. 
Equation  (3-3)  may  be  written  as: 


k=lj=k 
m  fi 


k=iM 
m  n 


k=lj=k  k=Ij=k 


(3-4) 
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Where: 


Remarks: 


1.  The  first  summation  term  of  equation  (3-4)  defines  a  symmetric  bilinear  form  or 
"sum  bilinear  form"  (sbf),  while  the  second  term  defines  a  skew-synunetric  or  the 
"difference  bilinear  form"  (dbf)  (also  known  as  the  oscillating  bilinear  form  [4]). 


2.  If  the  bilinear  form,  equation  (3-4),  is  solely  symmetrical  then: 


^jj< 


PjM 


W  _  ,['■]  _  ,['•] 


^j.k  -^kj’ 


41 


=  0 


And: 

(•.•)  =  Z  {^k«J  +  ) 

k=ij=k 

For  notational  consistency  and  to  facilitate  development  of  the  state-space  model 
utilizing  the  normalized  form  (to  be  considered  later),  we  shall  write  equation  (3-4)  as: 

X  llPk.jdkj(xkiij+XjU^)+  X  ll/^k,jd*k.l(x,,uj -XjU^)  (3-5) 

k=Ij=k  k=lj=k 


Where: 


Pk.j  = 


1  for  j^k 
1 


otherwise 


l  =  -\-j;  k-l,2,..,m;  j -k  + l,k  + 2,...,n 
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Remarks: 


1.  The  bilinear  form  of  equation  (3-5)  may  be  regarded  as  consisting  of  a  weighted 

sum  of  normalized  symmetric  elements  and  a  weighted  sum  of  normalized  skew- 
symmetric  elements.  Terms  |/?^.  ]|  will  be  referred  to 

as  the  normalized  symmetric  and  skew-symmetric  elements,  respectively,  of  the 
bilinear  form  while  coefficients  will  be  referred  to  as  the 

weights  on  these  elements. 

2.  Clearly,  k  =  m;j  =  «  =>  — m  +  l)  ^  ^  equation  (3-5),  the  number  of 


terms  in  the  sbf while  the  number  of  (non-zero)  terms  in  the 

rry  7  m{2n  -  m  - 1)  T.T  .7  .  ^  r  .7 

^ur  _  ;  ^ - i  Note  that  d  i,i  =  — - — -  =  0  for  j  =  k . 

Pkj 


dbf  =  l„,r„-ni  =  - 


3.  The  factor  has  been  introduced  in  equation  (3-5)  in  order  to  account  for 

"normalization"  of  the  bilinear  state-input  vector  (defined  below)  similar  to  the 
normalization  of  ihe  quadratic  state  vector  considered  in  section  2. 

3.1.  Definition 

The  "sum  hilinenr  state-input  vector"  (sbsiv)  will  be  defined  as  the  vector 

whose  elements  are  the  (un-weighted)  normalized  elements  of  the  symmetric  bilinear 

form  of  equation  (3-5),  that  is: 


=  1.2 . 


Where: 

PkJ  +  ’^j^k  ) 


3.2.  Definition 

The  "difference  bilinear  state-input  vector"  (dbsiv)  will  be  defined  as  the  vector 

whose  elements  are  the  (im-weighted)  normalized  elements  of  the  skew- 
symmetric  bilinear  form  of  equation  (3-5),  that  is: 
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Where: 

;r W* j (x^Uj  - XjUk ) 

/  =  (A:-/)j^M-|j  +  y;  j^k 


3.3.  Definition 

The  '‘bilinear  state-input  vector  function"  (bsivf)  (or  the  bilinear  forcing  function)  will  be 
defined  as  the  vector  whose  elements  are  the  (weighted)  normalized 

elements  of  the  bilinear  form  of  equation  (3-5),  that  is: 


. »} 


(3-8) 


Where: 

77;(x[^ImW)=  X  f,Pkjdk.i(xki*j  +  XjU^)+  X  'ZPk.jd*k.l{xkUj  -  XjUk) 
k=lj=k  k=Ij=k 

=  Z  Z  lLd*kj7rY^ 

k=lj-k  k=lj=k 


l^{k-l\n—^ 


+  J 


Remarks: 

1.  The  significance  of  the  symmetric  and  skew  symmetric  parts  will  become  clearer 
later  when  we  consider  the  basis  matrices  that  define  these  bsivs. 


3.4.  Example 

In  order  to  explore  certain  useful  properties  of  the  bilinear  state-input  vector  (bsiv)  or 
the  forcing  function  we  consider  equation  (3-5)  for  n=3,  m=2: 
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nf^  =  di,{x,Vi)  +  d,:,^{x,U:,  +  ,?/;)  + 1//, 5  dji-^ix^u,) 

d j,2  ^ —  ("^ 2  ^ 2^  1')  d  1^2  (  ^3 ^7 ) d j  2  (  ^3 ) 

/zl'^  =  d, , {x,Uj)+d^ , +  x,W; )+ C?, 3  -^{xjV,)+d2  J )  +  J, 5 ) 

■’■  V2  ’  v2  V/ 

+  dl  l  )+ ^ 

Tzl’^  —  d ^ j{x jU j')'¥ d ^ 2  {xjUi  +  ^2^ d 2^3  {xjUj )  +  <^3 ^ )■!' <^3 3  {x^Ui ) 

+  ^i?3  9  {^1^2  ~~  ■^9^7  )■*■  (~  X377y  )+  ^3,5  (“  ) 

These  equations  may  be  written  in  matrix  notation  as: 


TZj'^ 

-  [d;  ;  dy  9 

du  d,2  d„\^%,)+[dl2 

dlsV' 

''•(...) 

M'- 

^  =  \d  2,1  d2,2 

d2.3  d2,,  d2j\r^-\.,.)+[dl 

2  dls  dls\^'^'{.,) 

(3-10) 

^F' 

II 

_Ui 

‘tNj 

d3,3  d^^  ^(•’•)+ 

9  djj  ^55] 

k“'CJ 

or  alternatively  as: 

^wyfl„w)^[BMy9i^w,„w)+[BW*]2[-’r(5W,„i'i) 


(3-11) 


Where: 


[d'-]]={^/3,,.V-  =  7,2,5;/  =  U...,5}:  is  a  [ixj]  sum  bilinear  state-input  (sbsi)  coefficient 


matrix 


D^~^'\=^*,  :J  =  l,2,3;l  =  2,3,5}:  is  a  [ixi]  difference  bilinear  state-input  (dbsi) 
coefficient  matrix. 


Remarks: 


1.  Note  that£*^^(x‘^^,M^^^)  is  a  [5 x /]  s&siu: 
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s'V‘\u''')= 


(xjUj)  -^{xiUo+XjUi)  -^{xjUj)  (xjUj)  -^{xjU^) 


And  that  is  a  [5x7]  dhsiv: 

-^{xjU2-X2U,)  -^{-Xjll,)  -^{-XjU2) 


Also,  note  that  matrix  includes  only  the  non-zero  terms  of  {sf* , }. 

2.  In  the  above  example,  the  shsiv  may  be  regarded  as  a  vector  spanning 

the  space  whose  elements  are  the  terms  of  a  bilinear  polynomial  in 
arranged  in  lexicographic  order. 

3.  Similarly,  the  dbsiv  may  also  be  regarded  as  a  vector  but  spanning  the 

space  whose  elements  are  the  terms  of  a  bilinear  pol3momial  in 
arranged  in  lexicographic  order.  However,  (ignoring  those  terms  that  are  always 
identically  zero)  there  are  two  fewer  terms  involved  in  the  case  of  dbsiv. 

4.  The  inclusion  of  the  factor  — ^  caters  for  normalization  as  was  done  for  the  qsv  in 

42 

the  previous  section.  The  metric  property  of  the  sbsiv  was  considered  in  section 
(2.4.4)  where  it  was  shown  that  for  n=m: 


I  + 


JJ] 


.PJ 


(3-12) 


For  the  above  example  n=3,  m=2,  it  can  easily  be  verified  that: 


JP 


,1’] 


and 


JP\ 


(3-13) 

(3-14) 


Where: 


"  j=l 
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[/;]  = 


1  0 
0  1 
0  0 


5.  We  observe  that  a  bilinear  polynomial  (the  "sum"  bilinear  polynomial)  of  the  form: 


{xjUj)  +  ^{xjV2  +j:_,W;)  +  -^(x5M;)+  (x2M,)  + 
\2  V2 


may  also  be  written  as: 


(3-15) 


Where: 


k;]= 


1  0 
0  0 
0  0 


;  [fiA 


0  1 
1  0 
0  0 


■Afuhj; 


0  0 
0  0 
1  0 


/■  F2^2  .  ■ 


0  0 
0  1 
0  0 


0  0 
0  0 
0  1 


It  will  be  noted  that  if  both  n,  m=3  dien  the  matrices  ^  J  are  symmetrical  as  in  the 

case  of  the  quadratic  polynomial  considered  in  section  2.  Also,  these  matrices  are 
similar  to  the  basis  matrices  considered  for  the  case  of  quadratic  scalar  functions 
and  serve  the  same  purpose.  That  is,  they  enable  us  to  form  a  general  symmetric 

bilinear  scalar  function  in  systematic  way,  as  follows: 


I X  '£A;  [Fjj, Aj=a  scalar 


(3-16) 


6.  Similarly,  if  we  consider  a  bilinear  polynomial  (the  "difference"  bilinear 
polynomial)  of  the  form: 
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then,  this  may  also  be  written  as: 

=W'T|2;l:tei|y'') 


(3-17) 


Where: 


[0  7]  r  ^  ^  1 

o.'teb-t  0  »;fel=-4o  0  ;fe]=I"l 

0  0  ^^-10  0  -1 


It  will  be  noted  that  if  n,  m=3  then  the  matrices  [f*^]  would  be  skew-symmetrical. 

These  matrices  are  similar  to  the  basis  matrices  considered  earlier  for  the  case  of  the 
"sum"  bilinear  polynomial  and  serve  the  same  purpose.  That  is,  they  enable  us  to 

form  a  general  skew-symmetric  bilinear  scalar  function  p  in  systematic 

way,  as  follows: 

[.Mk-J  J  (3-18) 

/  =  (^-^)[«-|]  +  7V  k/l=M 

7.  As  discussed  in  section  2,  the  basis  matrices  may  be  used  to  establish  a  relationship 
between  an  sbsiv  and  the  corresponding  linear  state  vector  and 

linear  input  vector  as  follows: 

Now: 


{xjUj) 

1  (  \ 

K/IIl 

-j={XjU2  +  X2Uj  ) 
y]2 

hs] 

V 

(;c,W7) 

y^T  [F2.2] 

y^T  Ksi 

(3-19) 
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Xi  ^  r 
X2  X, 

yf2  42  L"-’ 
42 

0  x^ 

0  i 


mhiu  \z  :  : ![:: 

y'lf [f,,}  0  u2  0 

yffer  «  0  % 


Where: 


Xj  0 

Xj  Xj 

42  42 


=  43  0  ;and  U[i\  = 


VJ 

0  jc. 


Uj  0  0 

^2  U, 

42  42 

0  0  ^ 

42 

0  U2  0 

0  0^ 


u'li 


(3-20) 


8.  The  above  result  can  be  generalized  to  the  case  of  4^^  e  if" ,  4^^  €  if"’ .  Matrices 
e  and  C/|yj  e  ,  Ij  ^  -will  be  referred  to  as  the  sum 

bilinear  generator  matrices  (sbgms)  that  map  respectively  4’^  e  if"'  and 
e  if"  onto  if^^ ;  That  is: 


(3-21) 

It  can  be  shown,  that  provided  not  all  the  elements  of  w^^^and  4^^  are  identically 

r_[7]-|  r— r— [7]1 

zero,  then  there  exist  the  inverse  sbgms  (isbgms)  =  f7[yj 


such  that: 


Note  the  changed  super-  and  sub-scripts  in  the  inverse  generator  matrices. 


(3-22) 


Where: 


=1  x[/]  x\i]  I  xl/]  :  is  an  [mx/y]  inverse  sum  bilinear  generator 


matrix  for  4 
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r-wi 

f 

r-wf 

'r-wi 

[7^12]! 

i 

m 

V\i] 

} 

U[i] 

;  is  an  [wx/;]  inverse  sum  bilinear  generator 


matrix  for 


9.  For  our  particular  example,  it  can  be  shown  that: 


r 

1 

{2X]-  +X,"  +X2')  XjX2 

2([i] 

X[,] 

II  j 

J  ^ 

X,X2  (X;'  +  2x,'  +  X^-  ) 

T 

2 

U/  +  +  ^,"  )  -XjX-y 

Xli] 

Xij] 

\  J  J  3  }  12 

f  ^  2  ->  '>  \ 

[L  J 

J  ^ 

-XjXj  yxf  ^  Xn' -¥  x^  j 

Where: 

A  =  +  X2'  +  Xj'  \2xi'  +  2x,'  +  xf  )J 


And: 


-^[3] 


It.?  3\  x■,[2x^- +  Xs~j  Xjixj-  +2x-,~  +X2'j  I 

xj{xj-  +2X2-  +x/j  j 


{  2  i 

+l{2 

^2] 

-XjXiX^) 


j  Xj(2xj^^+x/)  x2(2x2^+X2^+xs^) 


V2 


10.  Similarly  it  can  be  shown  that: 


(3-23) 


—[2]' 


(  ’  ^ 

>  /  N 

7  M7" 

+“ir 

UjU2 

2 

V  ) 

K  2.  , 

M;M7 

0 


(2  ^ 
U]  0 

V  ^  ■ 


0 


\Uj  +-U2  j 
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T 

FttI-’]!] 

-1 


(  '> 

Uf  9 

I  ^  '  7 

(-WyW,) 

0 


(-W/W9) 


^  2\ 
9  W9 

+-T 

V  ^  y 

0 


2{uj'  +W9" ) 


Where: 

/!  =  {uj~  +1/9")" 


And: 


^[2] 


2Ui 


/  2 

W; 


+  M9 


(-MyM,') 

0 


2uy 

2uo^ 

vl 


{-UjUo^) 


2u7 


2u 


;(»/+?//) 


2^ 

9  W9 

w;  +-r 
^  y 


2m 


y(M;-  +M;-) 


r2 


r2 


(3-24) 


11.  Consideration,  analogous  to  those  discussed  in  remarks  7-10  above,  also  apply  for 
the  case  of  dbsiv,  the  results  are  quoted  below: 


2W(sl'l,.W)  = 


-^{xjU-> -XoUj) 
->12 


m 

pii'U''') 

m 

pi/] 

=r 

h/fe''’) 

= 

(,1'T 

Fj/] 

m 

^2,/] 

_£9_  ^ 

42  VJ 

-IL 


0 


Uj 

U2 


X 


[-’]* 

U] 


(3-25) 
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[Wl 

t' 

Fu. 

T 

y'Ti 

1 

Fo.s 

T 

Ho 

^2 

0 


JLl 

^2 


0 

0  0 


0 

JiL 

\r2 

U-, 

T2. 


X; 

^2 

- 

^\i] 

.^3. 

M 


Where: 


X 


[/] 


Xj 

Ho 

iij 

0 

r2 

r2 

X3 

~yr2 

0 

;  and 

^[/] 

= 

0 

0 

u, 

4~2 

0 

1 _ 

0 

0 

M, 

12.  The  above  result  can  be  generalized  to  the  case  of  e  R” ,  e  R”’ .  Matrices 


(3-26) 


X 


[/] 


e  and 


U 


[/] 


e  R^-^  will  be  referred  to  as  the  difference  bilinear 


generator  matrices  {dbgms)  that  map  respectively  ir’‘  e  R”‘  and  ^  ^  e  R” , 
1 2  =  — — — — ,  onto  R^- ,  That  is: 


X 


[^] 


,W- 


u 


I^] 


{2,-27) 


13.  Unlike  the  previous  case,  the  existence  of  an  inverse  generator  matrix  is  not 


jjW 


as 


guaranteed.  In  fact,  it  can  be  shown  that  an  inverse  dbgm  (idbgm) 
defined  by  the  relation: 

:  does  not  exist,  since  the  inversion  is  not 

However,  it  may  be  shown,  that  provided  that  none  of  the  elements  of  are 


f 

T 

r-u]*i 

1 

r 

'i 

I _ 

\ 

^[1] 

defined. 


identically  zero,  then  there  exist  an  inverse  dbgm  {isbgm) 


X 


[;]* 

[2] 


such  that: 


X 


[5] 


W 
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Note  the  changed  super-  and  sub-scripts  in  the  inverse  generator  matrices. 
Wliere: 

fr  .  .-iTr  .  ^.~\T 

:  is  an  [mx/,]  inverse  sum  bilinear 


1 

1 _ 

f 

T 

1^' 

= 

^[1] 

^[1] 

L  J 

1 

2 

J 

L  J 

generator  matrix  for  u 


[^] 


For  tire  current  example: 


^[2] 


’  (  2  2\ 
—  X2XJ  -X2\Xi  +X2  ) 


-XJX2X2 


XjXf 


■XJX2X2 


■2(3 


{x2^  +^3"). 


(3-28) 


Where: 

2I  2  2 

A  =  X2  \Xj  +X2  +^-5  j 

3.5.  Generalised  Bilinear  State-Input  Generator  Matrices 

For  a  state  vector  e  R"  and  an  input  vector  the  sbgms  were  defined  in 

section  2,  equations  (Bl-2),(Bl-3).  The  sbgms  and  dbgms  matrices  may  also  be 
constructed  utilizing  the  basis  matrices  J'  Following  the  procedure  given 

in  equations  (3-19),  (3-20)  of  the  example,  an  algorithm  for  constructing  the  basis 
matrices  [Fj  J  for  a  general  case,  may  be  defined  as  follows: 


for  j  =  k 


for  j  It  k 


1  for  r  =  s 
0  otherwise 
1 


for  r-k,s=jandr  =  j,s  =  k 


(3-29) 


V2 

0  otherwise 


j  =  l,2,...,m;k  =  j,j  +  l,...,n 

The  basis  matrices  for  a  general  n-state,  m-input  symmetric  bilinear  vector  system  are 
given  below: 
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Similarly,  the  dbsigm  may  be  constructed  using  the  basis  matrices  following  the 
procedure  given  in  equations  (3-25),  (3-26)  of  the  example,  an  algorithm  for 
constructing  these  basis  matrices  may  be  defined  as  foUows: 


for  r  =  k,s  =  j 

for  r  =  j,s  =  k 
otherwise 


(3-30) 


j  =  l,2,...,m;k  =  j,j  +  l,...,n;  andj^k 

The  basis  matrices  for  a  general  n-state,  m-input  skew  symmetric  bilinear  vector 
system  are  given  below: 
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The  dbgm  matrices  X  fjj  ,  U  fi]  are  given  in  Appendix-B,  equations  (Bl-4),  (Bl-5). 

3.6.  Bilinear  State-Input  Vector  Function 

Finally  before  we  leave  this  section,  we  give  a  state  space  representation  of  the  state- 
input  vector  (or  a  forcing)  function.  Using  relationships  developed  in  this  section,  a 
bilinear  n-state,  m-input  vector  function  may  be  written  as: 

Mv  -  «''!)=  .l'>)  (3-31) 


Or  alternatively,  in  terms  of  the  hgms,  as: 


(3-32) 
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4.  Forced  Quadratic-Bilinear  Dynamical  System 

We  shall  now  consider  the  algebraic  structure  of  a  forced-quad ratic/bdinear  system 
where  the  RHS  of  the  differential  equations  consists  of  linear  and  quadratic  terms 
identical  to  those  in  the  free-QDS  as  well  as  additional  linear  input  terms  and  bilinear 
terms  involving  inputs  (forcing  functions  or  control  inputs)  and  system  states. 


4.1  Definition 


A  dynamical  system  will  be  classed  as  a  forced-quadratic/ bilinear  dynamical  system 
(forced-QBDS)  if  the  right  hand  side  (RHS)  of  the  differential  equations  describing  die 
time  evolution  of  states  consists  of  (see  sections  2  and  3): 

a.  a  linear  term  in  system  state  vector: 

b.  a  quadratic  term  in  system  state  vector: 


c. 

d. 


a  linear  term  in  input  vector: 
a  bilinear  term  in  state-input  vector: 


That  is,  the  forced-QBDS  may  be  written  as: 


+ [fii-’iy-’i + kwy'i +awy'i,„w) 
= U'ly'i + + [cwy'i + [Di-’ikwy'i,„w)+ 


(4-1) 


Where: 


=  =  7,2,..., «}  :  is  an  [w  x  «]  linear  state  coefficient  matrix 

=  W . ^ 


1  . 

nin  +  7) 

> :  is  an 

n  X— - - 

2 

quadratic  state  coefficient 


matrix 


=  {cj  ^ ;/  =  ],2,...,n;k  =  l,2,...,m  <  «}:  is  an  [«  x/m]  input  coefficient  matrix 


djj  :  j  =  l,2,3,...,n;l  =  1,2,3,..., 


m{2n  -  m  + 1) 


IS  an 


wx 


m{2n-m  +  ]) 


sum 


bilinear  state-input  (sbsi)  coefficient  matrix 


1 

1 

£,l-J 

=  \d jj :  j  =  1,2,3, ...,n;l -2,3,5,...,  ^  h 

is 

an 

«x - 

2 

difference  bilinear  state-input  {dbsi)  coefficient  matrix 
=  7,2,...,«}:  is  an  [«  x7]  linear  state  vector 

n{r}  +  l) 


=  1,2,..., 


n{n  +  l)\  . 

'  ■  ^ :  isan 


-x7 


quadratic  state  vector 


=  l,2,...,m  <  «}:  is  an  [/m  x  7]  linear  input  vector 
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n}-  =  \n]~^;i  =  is  an  x  /]  bilinear  state-input  vector  function 


(2/7  -m  +  7)1  _  . 


IS  an 


m{2n  -  m  + 1) 


x7 


sum  bilinear  state- 


input  {shsi)  vector 

. . 

bilinear  state-input  (sbsi)  vector 


IS  an 


m{2n-m-l) 


x7 


difference 


Alternatively,  utilising  the  qgm  and  bgm  defined  in  sections  2  and  3,  the  Forced-QBDS 
may  also  be  written  as: 


X 


[/] 


+ 


D 


Or  equivalently  as: 


(4-2a) 


(4-2b) 


Tlie  forcing  function  or  the  input  vector  function  Ui=  (wy  ;i  =  7, 2,. ..,77}  is  given  by: 


(4-3a) 

(4-3b) 


(4-3c) 


Where: 


[/]  J’ 

r— i-’]i 

— I 

I 

_ I 

'-tM' 

’77W*’ 

rwj’ 

■  ^[1] 

:  are  the  qgm  and  bgms  defined  earlier  in  sections  2 


and  3. 


4.2.  Example: 

For  n=3  and  m=2,  the  forcing  function  or  the  input  term  Wy  = 
may  be  written  as: 
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+ ‘^1,2^2  +di  jX,Uj  +djj-^{xiU2+X2Ui)+di3  -^{x JU i)+d J^X 2^ 2 

+d  ]  j  -^{x2U2)+d  ]  2  ~^2^*])'^1,3 

Up^  =  C2jUj  +C2  7«2  +d2jXjUj  +d2  7  ~^{^I^2  +  ^ 2^ l)'^ 2,3  ~1=  3’’'' l)'^ 2J^ 2^h 

'  V2  '  \2 

+d2^S  -^{x3U2)  +  d2j-^{xjll2-X2U])+d2.3^i~^3^<l)+^2.5  -^(-^3^2) 
-  ^3,1^1  '^^3,2^2  '^^3,1^1^I  '^d32~^{xiU2  X2U  i)+d  3J -^{x  3U  i)+d  j^X  2U2 

^  ’  ’”  ’  ’  V-2  \2 

+d3  3-^[x3V->)  +  dl^^{xiU2-X2U,)+dl3-^{-X3Ui)+dl3^{-X3V2) 

’  V2  V2  ^2  yl2 


Remarks:  Compare  the  above  result  with  equation  (3-11). 


5.  Quadratic-Bilinear  System  Output  Model 

Utilising  the  definitions  given  in  previous  sections,  dynamical  system  output  or 
measurement  models  may  be  derived.  Following  on  from  the  results  of  previous 
sections,  the  system  output  model  may  be  written  as: 

jW^[„W],Mq^W],W  +  „wy'l.,W)+[iW]vH  (5-ia) 


(5-lb) 


Or  using  qgm  and  bgms  as: 

.l'>={[/,l'l]+kW][A-JSlh['l+|[iW]+[MWF-« 


+ 


X 


[/] 


M 


[2]* 


X 


[/] 


.w 


^[1] 


[2] 


(5-2a) 


(5-2b) 


Where: 


[H^’^\={hij;i  =  I,2,...,r;j  =  ],2 . n}  :  is  an  [r  x«]  linear  state  output  coefficient  matrix 


is  an 


n{n  +  i) 


quadratic  state  output 


coefficient  matrix 

{/,  t,-/  =  l,2,...,r;k  =  l,2,...,s  <  «}:  is  an  [r  X5]  output  disturbance  coefficient  matrix 
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s{2n  -  s  +  ]) 


{2n-s-])\  .  r 

IS  an  rx 


s{2n  -  s  - 1) 


\  ,f’ll  f  ,-1  s{2n-s  +  l)]  .  r  s(2n-s  +  l) 

[M'-'\=  =  l,2,...r; p  =  1,2,...,-^ — - ^4:  is  an  rx - - - 

state-disturbcince  (sbsd)  coefficient  matrix 

‘  f  *  .  ,  „  ,  „  s{2t7-s-l)]  .  r  s{2n-s-. 

=\mip;i  =  l,2,...r;p  =  l,2 . — - - '4:  is  an  rx - - - 

bilinear  state-disturbance  (dbsd)  coefficient  matrix.  Note  that  matrices 

are  similar  to  I  of  the  previous  sections  3  and  4. 

=  l,2,...,r  :  is  the  [r  x  /]  output  (measurement)  vector 
=  l,2,...,s  :  is  the  [j  x  /]  output  disturbance  vector 


sum  bilinear 


difference 


A"|)]  ;  :  are  respectively 


s{2n-s  +  l)  ^ 
2  ' 


X5  and 


s{2n  -s  -]) 


X  5  bgms 


containing  elements  of  4^^  similar  to  those  have  been  defined  in  sections  2  and  3. 


-  7  -  7 


V  ;  V  [)j  :  are  respectively 


s{2n  -  s  + 1) 


xn  and 


s{2n-s-l)  ^ 


xn  bgJHs  containing 


elements  of  similar  to  those  defined  in  sections  2  and  3. 


A  block  diagrams  for  the  forced-QBDS  with  output  model  is  given  in  Figures  5.1. 
Figure  5.2  gives  a  block  diagram  where  control  input  and  disturbance  vectors  are 
linear. 

5.1.  An  Example  of  a  Second  Order  Matrix  Riccati  Differential 
Equation: 

We  consider  the  matrix  Riccati  differential  equation  of  the  type  commonly  encountered 
in  synthesis  of  an  optimal  controller  or  a  filter  for  a  linear  dynamical  system  [9,10].  The 
matrix  Ricatti  differential  equation  resulting  from  the  estimation  (filtering)  problem  for 
a  second  order  system  is  as  follows: 

4  [^-l = 1^1M+ UY  W-  W[i]  W+ to]  (M) 


Where: 


:  is  a  [2x2]  linear  state  coefficient  matrix 

^2,1  ^2,2 


<11  92  , 


92  93 


:  is  a  [2x2]  (symmetric)  system  noise  covariance  matrix 


Pi  P2  , 


P2  P3 


:  is  a  [2x2]  (symmetric)  Riccati  matrix 
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^1.1  ^],2 
Sj2  ^2,2 


:  is  a  [2x2]  (symmetric)  matrix 


Expanding  equation  (3-33)  in  its  elemental  form  gives  us: 
^Pj  =  2(ajjpj+aj2P2)- \^i,iPi'  +  2s  12P1P2  +S2^2P2'\+<ii 


dt 


^P2=a2iPi-^{aji+a2i)P2'^‘^2.lP3  ~  ^1.2P2^  +^l.lPlP2  +^l,2PlP3  ^  2.2P2P3  + 2 


dt 

^P2-2{a:,jp-,-^a-,-,p3)-\sj,po'  +2si  -f-5, ,/?/]+ 

dt 


(5-4) 


Equation  (3-34)  may  be  written  in  matrix  form  as: 


Pi 

2^1,1 

2^1  y 

0 

Pi 

P2 

= 

^2,1 

{<3 1 1  +  02,2 

)  <^2.1 

P2 

P3, 

0 

2‘^2,1 

1 

<N 

^P3^ 

which  is  of  the  form: 

=  \aI“]x"’  +  [s® 


Pi 


3i,i 

3i,2 

0 

3  2,2 

0 

0  ' 

P1P2 

^1 

0 

^1.1 

^1,2 

^1.2 

^2,2 

0 

P1P3 

2 

+ 

^2 

0 

0 

0 

^1,1 

Si  2 

^2,2 

P2 

.^3_ 

P2P3 

■? 

.  P3'  \ 

(5-5) 

(5-6) 
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6.  Extension  to  a  Second  Order  Approximation  of 
Analytical  Functions 

One  obvious  application  of  the  quadratic/ bilinear  algebraic  structure  discussed  in  the 
previous  sections  is  in  the  second  order  approximation  of  analytical  functions.  Non¬ 
linear  dynamical  systems  witli  RHS  that  consists  of  analytical  functions  can  be 
conveniently  set  up  as  a  set  of  quadratic/ bilinear  differential  equations  [5].  In  order  to 
demonstrate  this  we  consider  the  following  non-linear  dynamical  system  of  the  form: 

^  Xj  =  fi  (x,  m);  i  =  1,2,  ...,n  (6-1) 

at 

with  the  output  model  given  by: 

Zj  =  hj{x,v);  j  =  l,2,...,r  (6-2) 


Where: 

X  =  {xjii  - 1,2,.. .,t7}:  is  an[«  x7]  linear  state  vector 
u  =  {uj  ;i  =  l,2,...,m  <  w}:  is  an  [w  x  l\  linear  input  vector 
y  =  {v,  ;/  =  l,2,...,s]-.  is  an  [s  X  /Joutput  disturbance  vector 
z  =  {zj-.i  =  l,2,...,r]\  is  an  [r  x l\  output  (measurement)  vector 


Assuming  that  functions  {/)(.,.)}  and  {/?,(.,.)}  satisfy  the  usual  conditions  for  the 
existence  of  the  Taylor's  series,  it  follows  that: 


dx 

du 


Sx  =  Y, 

)  j=l 

m 

I  J=1 


J 


du, 


dx 


dx 


A 


dx 


du 


du 


^5u 

—  V 


I 

,  k=i  ^Xk  \^j=i 


—  Sx, 

ydXj  j  j 


,  k=l  j=I 


z 

I  k=I  ^^k  J=I 


(  A  ^ 


du  j 

^dUj)  ^ 


,  k=l  j=I 


'  a-/  ' 

dXkdXj^ 

dx.duj^ 


dx 


\  ^x 


^du 

ydu 


r  „  \ 


du=Tz^  y 

,  k=i  ^^k  y  j=i 


5/ 


J 


\dxj 


,  k=I  j=l 


”  ^  aV/  ^ 


ydu^dXj  J 


y 


\duj 


(  d'ft  '' 


'.=ZZ 

k=I  j=I 


^dUkdUj 


dXk  dxj 


dxj^duj 


dxjdUk 


du^  du  J 
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Since: 

n  m 

sz 

k=l  J=I 


'  s-f,  ''' 

ydx^duj  j 


III  /I 


k=i  j 


Sx  jSiii^ 


the  second  order  Taylor's  expansion  gives: 


j=r 


"  ^  dfi  ' 


^'/  +  E 

7=/' 


^  dfi  ^ 


y^dUj  J 


d-f,  ^ 

dxjfdxj 


dxj^&c  j 


Ao 


+  ZZ 

k=lj=l^ 


dxi^du  I 


^T^,/+-rZZ 


a-y;.  ^ 


du^dii 


k^lj=K 

dukSufi  i  =  l,2,...,n 


J  J 


(6-3) 


Similarly: 


dXi 

n  s  I  3-'/;, 


&,+x 


^  i"  a//, 


7  -^3,  ^ 


av. 


inn 


+  ZZ 

k=lj=l\ 


[^dxkdvj^ 


2k=ij^]\ 


K''-JJ\0  ~k=Ij=l^ 
-^3,  ^ 


Z'  32.  ^ 

- —  SXk&Cj 

^BXkdXj)^ 


a-a, 

v^VA-dv,-, 


dv^dv.;  i  =  l,2,...,r 


(6-4) 


Furthermore  since: 


ZZ 

k=l  j=l 


'  a-/,  ^ 

^ax^aw^.  y 


,  A=i  7=A 


\{  a-’/  1 

dxiSiii  -t- 

f  a-7.  1 

&(  i&t^  1 

[lax.awj 

A  J 

0 

dXjdUk  J 

. '  J 

^k.j=S 


for  k^j 


1  otherwise 


It  follows  that  equations  (6-3)  and  (6-4)  may  be  written  as: 


f{x,u)  =  f(x,u\^+Y. 

j=r 


"  ^  d/f 


^  Y  . 

Ao 


^j  +  T 

j=t 


m  f  ay;.  M 


\^^JAo 


&lj  +  ^ 

k=lj=k 


^^k^^JAo 


Sx^Sx  j 


+ 1  'L<^k.j 

k=Ij=k 
i  =  l,2,...,n 

And: 


^  a-y;- 


dx^du 


JAo 


Sx^Suj  + 


^  a-y;-  ^ 

^dXjdUk  y|^ 


dXjSUk  [  +  ££«^y,A 
k=iM 


J!a 

du^du 


J  J 


SukSuj 


(6-5) 
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^  a/?. 


dh^ 


/?,.  (x,  v)  =  hf  (x,  v)|^  +Z  T-^  +  Z  ^  ‘^y  +  Z  Z«^A',y 


—  Sxj^Sxj 


+  ZZ«A,y  ^ 

k=ij=k  l^l^ox^dw 
i  =  l,2,...,r 


dKf-dVj  + 


dxjduj^ 


/■Si'  3-/7. 

+7ZZo^A,/ 

,  ■'  ^A=/y=A  [dv^dVj 


dvj^dvj  (6-6) 


6.1.  System  Dynamical  Model 

Making  the  following  substitutions: 

,  j  =  l,2,...,n;  i  =  l,2,...,n  ■,  Cjj  =  ;  j  =  l,2,...,m;  i  =  l,2,...,n 


1  f  ']  7  7-1  nin  +  l)  .  1/77/  ^ 

bij  ^  x'  ;  1  =  1,2,...,-^ — ;  /  =  7,2,...,/7; /  =  {A:-7)  ///-- 


dx/^dxj 


r  1  r  r)~  f  ^ 

j  =  l2,...,m;  k  =  j,j  +  ],...,n 


dx  jduj^ 


duj^duj 


Using  the  relationships  of  equations  (6-3),  (6-1),  The  second  order  perturbation  model 
(differential  equation)  for  system  dynamics  may  be  written  as: 

J  77  77)  77  77  77  777  C  r,  f  1  1 

—  =Z«/,y^y  +  Zc,jA  +  Z  Z«y.A*7.A‘^A^y  +  Z  Z«^A,y  tl'y^AA  +4a^7^'a} 


k=ij=k 


Where: 


jU  =  {k-lin--\  +  j;  T]  =  ik-jim--  +J 


The  first  two  terms  on  the  RHS  of  equation  (6-7)  are  linear  in  {^}  and  {^} 
respectively,  while  the  third  and  fifth  terms  are  quadratic  involving  and 

respectively.  The  fourth  term  is  the  bilinear  state-mput  term  in  \a,a}.  Utilising  the 
results  of  previous  sections,  equation  (6-7)  (after  some  algebraic  manipulations)  may  be 
written,  in  state-space  form,  as: 
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dt 


Or  alternatively  as; 


dt 

+  fz)i-TVl-’l*|&l'i,&l'l)+[cl-’l]&M 


6.2.  System  Output  (Measurement)  Model 

Making  tlie  following  substitutions: 


hn=  ;  j  =  l,2,...,n;  i  =  l,2 


,2,...,r-,hj=\^\;  j^l,2,...,s;  i  =  l,2,...,i 


(  7)-h  '' 

4'1=  -''  ;  j  =  l2 s;  k  =  j,j  +  l,...,n 

44  ^  J  ’  j^I<2,...,s;  k  =  j,j  +  l,...,n 

«//=  ;l^l,2,...,^^^^;i  =  l,2,...,r-,  I  =  ik-lis-^]  +  j 

■  [dv^dvj )  2  V  2J 

Using  the  relationships  established  in  equations  (6-3)  and  (6-1),  the  second  order 
perturbation  model  (measurement  model)  for  the  output  may  be  written  as: 

=thjSKj  +  +  X  t^j,kji,f^&CkSKj  +  X  t^kAkJ^k^J  +fjl^j^k  } 


k=lj=k 


+  Z  llfXi.k^i.n^k^j;  i  = 

k==ij=k  ■ 


(6-10) 


Where: 

M  =  {k-lin-^  +  j;  ri  =  {k-lis-A  + j 


The  first  two  terms  on  the  RHS  of  equation  (6-10)  are  linear  in  {&}  and  \Sv\ 
respectively,  while  the  third  and  fifth  terms  are  quadratic  involving 
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respectively.  The  fourth  term  is  the  bilinear  state-input  term  in  Utilising  the 

results  of  previous  sections,  equation  (6-10)  (after  some  algebraic  manipulations)  may 
be  written,  in  state-space  form,  as: 

(6-11) 

Or  alternatively  as: 

&|'1 = Iff  H-ljWj&W + 


Figure  6.1  Block  Diagram  of  the  Forced-QBDS  of  equations  (6-9),  (6-12) 
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7.  Conclusion 

In  tills  report,  formal  definitions  of  free  and  forced  quadratic/ bilinear  dynamical 
systems  were  given  and  the  algebraic  structure  of  this  class  of  systems  was  explored 
witli  a  view  to  setting  up  a  systematic  approach  for  deriving  state  dynamics  and 
measurement  models.  This  class  of  dynamical  systems  is  characterized  by  a  set  of  first 
order  differential  equations  with  the  RHS  that  contains  linear  and  quadratic  terms  in 
system  states  as  well  as  bilinear  terms  involving  state  and  input  (or  control)  variables. 
An  example  of  the  quadratic  state-space  model  of  the  airframe  aerodynamics  is  given 
in  Appendix-A.  Properties  of  the  quadratic  and  bilinear  vectors  were  investigated  and 
relationships  between  these  and  linear  vectors  established.  Systematic  procedure  for 
constructing  quadratic  state  and  the  bilinear  state-input  vectors  was  derived.  The 
quadratic  algebraic  structure  was  applied  to  the  formulation  of  a  state-space  model  for 
a  second  order  approximation  of  a  general  (analytic)  non-linear  system. 

The  concept  of  quadratic  and  bilinear  generator  matrices  was  introduced  that  allows 
linear  (state  and  state-input)  vectors  to  be  mapped  onto  the  quadratic  and  bilinear 
vectors.  Properties  of  these  generator  matrices  were  explored  and  it  was  shown  that 
inverse  generator  matrices  may  be  defined  (in  all  but  one  case)  that  allows  quadratic 
and  bilinear  vectors  to  be  mapped  onto  linear  state  and  linear  input  vectors.  State- 
space  representation  of  the  system  output  or  measurement  model  was  also  derived  for 
the  case  where  this  contains  linear  and  quadratic  terms  in  states  and  bilinear  terms  in 
state-disturbance  (noise)  variables. 

The  state-space  structure  of  quadratic/ bilinear  dynamical  systems  considered  in  this 
report  should  facilitate  analysis  of  this  class  of  non-linear  systems  and  possibly  lead  to 
general  synthesis  techniques  or  extension  of  linearised  techniques  for  applications  to 
these  problems. 
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Appendix-A 


Euler's  Equations  of  Motion 

The  six  equations  of  motion  for  a  body  with  six  degrees  of  freedom  may  be  written  as 
[Ref.  5]: 

m(v-\-wq-vr)  =  X  -^T  + 

m(v  +  ur-wp)  =  Y  +  gym  (Al-1) 

m(w  -  ug  +  vp )  =  Z  +  g^m 

hxp  -  ( lyy  -  4- kr  +  Jyz ( ^  - 1 zx( P^i  +  ^ + 1 ^,(rp  -  q)  ^  L 

Iy^,q  -  ( /„  -  )rp  +  (p-  -r^  )- I^.(qr  +  p)  +  (pq-r)  =  M  ( Al-2) 

Jzzr -(Jxx-  lyy )P^  +  ^X)’(^^  -P^)- Jyz (fP  +  ^)  +  4 (qr-p)  =  N 


Separating  the  derivative  terms  and  after  some  algebraic  manipulation,  equations  (Al- 
1),  (Al-2)  may  be  written  in  a  vector  form  as: 
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q' 

qr 

■) 

r 

(Al-3) 


(Al-4) 


Note  that  the  states  (u,v,w,p,q,r)  appear  as  quadratic  terms  in  equations  (A  1-3),  (Al-4). 
Where: 
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yz^  xy  ^  yy^  zx  . 

+  44)  (44-4'A 

A  =  [ 

)jy 

,1- 

L4'-2444) 

Combining  equations  (Al-3)  and  (Al-4),  we  obtain  the  full  6*  order  rigid  body 
dynamics  state  equations  as: 


d 

i/ 

■[c] 

[0] 

_  j _ 

■[/] 

I  [0]  ■ 

1 

g 

dt 

±2 

1  I^J-'lsll 

±2 

+  ■ 

[0] 

j  |.r'l 

\y 

“T 

0 

Where: 

'0  0 


[C]  = 

=  [m  V 


0  10-1 
0-1001  0 
10-1000 

Y-,xf  =  [p  q  rf 
=  \uq  ur  vp  vr  wp  wqY 1^2^  =  \p^  pq  pr  q' 

^l^  =  [x  +  f  Y  Zj-,^?  =  [L  M  =  gy  gj 


(Al-5) 


Equation  (Al-S)  may  be  written  in  a  compact  form  as: 

dt~ 


Where: 


:  is  the  6x12  quadratic  state  coefficient  matrix 
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12x1  quadratic-state  vector 

I  =[x  +  7’  Y  ZjZ  M  TV^]^  :is  6x1  a  vector  function  of  control 


inputs,  forces  and  moments 

=  k i  of  =  k.v  Sy  Sz  o\o  of:  is  the  6x1  gravity  (or  disturbance)  vector 


45 


DSTO-TR-1497 


Appendix-B 
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;  is  the 


m 


{2n  -m-l) 


X  n 


m-input  difference  bilinear  generator  matrix  (dbgm). 
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